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ABSTRACT 

We search for static solitons stabilized by heavy fermions in a 3+1 dimensional Yukawa 
model. We compute the renormalized energy functional, including the exact one-loop 
quantum corrections, and perform a variational search for configurations that minimize 
the energy for a fixed fermion number. We compute the quantum corrections using a 
phase shift parameterization, in which we renormalize by identifying orders of the Born 
series with corresponding Feynman diagrams. For higher-order terms in the Born series, 
we develop a simplified calculational method. When applicable, we use the derivative 
expansion to check our results. We observe marginally bound configurations at large 
Yukawa coupling, and discuss their interpretation as soliton solutions subject to general 
limitations of the model. 
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1. Introduction and Motivation 


Chiral gauge theories, such as the electroweak Standard Model, present a challenge to 
conventional notions of decoupling. By increasing its Yukawa coupling, one can make a 
fermion so heavy that it should be irrelevant to low-energy physics. On the other hand, it 
cannot simply disappear from the theory, since then anomalies would no longer cancel. It 
is known that decoupling a chiral fermion leaves behind a Wess-Zumino-Witten functional 
of the Higgs and gauge helds, which keeps the path integral gauge invariant |^. However, 
to cancel Witten’s non-perturbative SU{2) anomaly |^, an even number of chiral doublets 
must be present in the low-energy theory. After decoupling a fundamental fermion, the 
presence of a fermionic soliton in the low-energy theory would ensure that gauge invariance 
is maintained at the level of the states [|l|. 

There is a natural mechanism for realizing this scenario. A twisted conhguration of the 
Higgs held will cause one fermion level to become very tightly bound, or even to cross zero 
energy Q . In the former case, the level can be hlled at very little cost in energy, while in the 
latter case the background held itself carries fermion number. This energy must be added 
to the classical energy required to twist the Higgs held, and then compared to the mass of 
a fermion in an unperturbed background. However, one must also include the contribution 
of the shift in the zero-point energies of all the fermion modes, since it is of the same order 
in h as the energy of the hlled level. If the total energy of a fermion number one twisted 
conhguration is below the mass of a free fermion, then the conhguration is stable. 

In this paper, we explore this phenomenon in a simplihed version of the electroweak 
sector of the Standard Model. We consider a Higgs doublet chirally coupled to a single 
heavy fermion doublet. We assume the fermions in the doublet have equal masses, and 
consider a hedgehog conhguration for the Higgs held. We ignore the SU(2)i-gauge helds. As 
we will argue in our conclusions, this may be an important omission, and work is currently 
underway to extend the calculation to include that case. 

Previous work |Q showed that quantum-stabilized chiral solitons do exist in the one¬ 
dimensional analog of the model we are considering. Here, as in that work, we will consider 
only the quantum correction from the fermion loop, which we expect to be the most impor¬ 
tant ehect. Neglecting bosonic loops is rigorously justihed in the large N limit, where N is 
the number of independent fermion species, but here we couple to only a single doublet and 
take N = 1. 

Our methods allow us to evaluate the one-loop fermion contribution to the energy exactly, 
maintaining a hxed renormalization scheme while exploring diherent Higgs backgrounds. 
There have been earlier attempts to compute the one fermion loop energy for a three- 
dimensional chiral background within various approximation schemes. Examples include 
discretization methods Bi, which may not be rigorously valid in the continuum limit and 
require cutoffs, making the renormalization obscure. Other approaches use expansions that 
are valid for slowly |0] or rapidly |^, |[ varying background helds. Also truncations 0 in 
the potential generated by the background held, the heat kernel expansion [M, and subsets 
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of field configurations with less severe ultraviolet divergences have been studied. While 
these methods may be appropriate for specific applications and useful within certain regions 
in configuration or parameter space, they do not allow one to explore a full range of ansdtze 
for the Higgs background or to make dehnitive statements about the existence of a soliton. 

We include in this Introduction a brief review of our method for computing the contri¬ 
bution to the energy from vacuum polarization induced by the background field. For a fuller 
discussion of the method see Ref. B.II1. 

The vacuum polarization energy is given formally by a sum over bound states and an 
integral over continuum energies in the background given by 




( 1 ) 


where {ej} are the bound state energies, uj{k) = yW+m? is the energy of a scattering state, 
and dn/dk is the change the continuum density of states due to the background held, dn/dk 
can be calculated from the derivative of the phase shifts with respect to fell. 


— =-trlnS'(A;) = —77 D{i)6{k,£) 

dk 2m ^ ’ vr dfc V 


( 2 ) 


where we have expanded the S matrix in partial waves labeled by £, representing quantum 
numbers like angular momentum and parity. d{k,£) gives the phase shift and D{£) the 
degeneracy. Since we will consider cases where the spectrum is asymmetric in energy, we 
dehne 6{k,l) to be the sum of the phase shifts over both signs of the energy: 6{k,l) = 
5+(/c, 1) + 5-(A:, 1), where ± specihes co = ±\/W+rn?. 

Before we can do the k integral, we must deal with potential divergences. We regulate the 
theory by dimensional regularization, that is, we analytically continue the entire theory to 
n dimensions, where the integrals converge. For integer n, the expected divergences emerge 
from the high momentum behavior of the phase shift integral. At high momentum, the 
Born series becomes a good approximation to the phase shift, so by subtracting successive 
terms in the Born series from 5{k, /), we can remove terms from the ^-integration that would 
diverge at physical values of n. The expansion in the Born series can then be unambiguously 
identified with the expansion of the effective energy in terms of Feynman diagrams with 
insertions of the background held 0 |l^ , [T^ . In noninteger dimensions, the contributions 
of both the Born terms and the Feynman diagrams to the vacuum polarization energy are 
hnite and unambiguous analytic functions of n. Therefore, when we subtract a term in the 
Born expansion and add back the equivalent Feynman diagram, we can be certain that we 
are not introducing hnite ambiguities into the computation. The subtracted integration over 
the density of states is then hnite and we can take the limit to integer n without difficulty. 
Finally, we introduce the contributions from the counterterms, which have been computed 
using standard renormalization conditions in the perturbative sector of the model. As usual, 
the potentially divergent pieces of the Feynman diagram are canceled by the counterterm 
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contributions. In all, the renormalized vacuum polarization energy is given by 

_ POO rjh _ rj _ / ^ma x \ \ 

-^max 

+ y r('^)(A) + rct(A) (3) 

N=1 

for bosons and fermions respectively, where the A^**^-order Born approximant to the phase 
shifts is denoted by 5^^\k,l). A^max is the number of Born subtractions required to render 
the k integration hnite. The compensating Feynman diagrams are denoted by r*^^i(A), 
and rct(A) represents the contribution of the counterterms. Both are cutoff-dependent, but 
as usual, renormalization conditions will determine an unambiguous, hnite result. We are 
left with two hnite and numerically tractable objects, the momentum integral and the sum 
Z]iv=T + ^ct- As a result, no explicit cutoh needs to be introduced in the numerical 
computation. 

This paper is organized as follows. In Section 2 we briehy review the Higgs sector 
of the standard model and outline its connection to the Yukawa model. In Section 3 we 
discuss the renormalization of the fermion loop and describe our calculation of the associated 
contribution to the energy. Section 4 contains the numerical analysis. We summarize and 
provide an outlook on future studies in Section 5. Technical details are given in three 
Appendices. In Appendix A we explore the Dirac equation and its scattering solutions for 
a chiral background held. In Appendix B we describe and numerically verify a simplihed 
treatment of contributions to the renormalized vacuum polarization energy that are higher 
order in the background held. In Appendix C we derive results in the derivative expansion, 
which we use to check our results in the case of slowly varying background helds. 


2. The Model 


The model we consider consists of the Higgs sector of the Standard Model coupled to a 
fermion doublet in 3 -|- 1 dimensions. Our goal is to explore the possibility that within this 
model there is a non-trivial Higgs held conhguration with nonzero fermion number whose 
energy is less than that of a state with the same quantum numbers built on top of the 
perturbative vacuum. The fermions get their masses through their Yukawa coupling to the 
Higgs. Our model dihers in two essential ways from the Standard Model: we omit gauge 
helds and our fermions have equal masses. At the end of the paper we discuss the possible 
sensitivity of our results to the omission of gauge helds. 

We write the Higgs sector of the Standard Model in terms of a 2 x 2 matrix-valued Higgs 


held 

V <^+ ) 

where {ipQ, </?+) is the usual doublet. The Higgs Lagrangian is 


(4) 


£H = -tr 




-Y(4>) 


(5) 
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where 


A 


V{ip) = — (^tr 


We take the vacuum expectation value to be 


($) = u 


1 0 
0 1 


and note that the Higgs particle has mass mn = \fXv. 

The coupling to the fermion doublet q = {t, b) is given by 


= gqi^Qn + gqn^'^qL 


( 6 ) 

(7) 

( 8 ) 


which results in mass m = gv for both t and b. It is also convenient to rewrite <h in terms of 
four real (dimensionless) helds s and p as 


^ = V {s + if ■ p) 


which gives 


/Ihf = mq{s + i'f^f ■ p) q 
With these dehnitions, the classical energy is 


^ci[*h] = (disdiS + dip- diP + - 1)^^ . 


(9) 

( 10 ) 


( 11 ) 


3. The Fermion Loop 

In this section we discuss the contribution of the fermion vacuum to the total energy. 
This contribution arises because the fermionic vacuum is polarized by the Higgs background. 
In order to compute this contribution we hrst have to outline the renormalization process 
in the perturbative sector of the model. The divergences of our model can be canceled by 
counterterms of the form 

Cct = atr — 6tr — ctr (12) 

where a, b, and c are cutoff-dependent constants. The Yukawa coupling g, and consequently 
the fermion mass m, are not renormalized at this order. 

In terms of the shifted Higgs held h = s — v, our renormalization conditions are that 
the vacuum expectation value of h vanishes, and that the fermion loop changes neither the 
position ruH nor the residue of the pole in the two-point function for h. In order to hx the 
counterterms, it is therefore sufficient to expand^ 

‘5eff[h] = —i Tr In {i^ — g{y + h)} (13) 

^Here and in what follows tr refers to sums over discrete labels while Tr includes the space-time integration. 
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up to quadratic order in h and combine the result with / d^xCct- In dimensional regulariza¬ 
tion we obtain 


a = 



— 6 / dxx(l — x) In 
Jo 


1 — x(l — x) 


im? 


mi, 3 ri , , 

-—r-/ dx In 

Am^ 2 Jo 


1 — a;(l — x) 


m 


H 




(14) 


where d = 4 — 2e, /i is the scale required to keep g dimensionless. 

Having set up the model in the perturbative sector, we now turn to non-trivial held 
conhgurations. We restrict our attention to the spherical ansatz for the Higgs held, 


^{x) = V [s{r) + if ■ Xp{r)] (15) 

with r = \/W^. With the standard form of the Dirac matrices, the corresponding Dirac 
operator becomes 

ms{r) —id-'V + mif-xp{r) 

—id-'V — mif-xp{r) —ms{r) 

and the fermion held obeys the time-independent Dirac equation. 




hnA/ = ooAl . 


(17) 


Note that the energy eigenvalue oo can assume both positive and negative values. In general 
the spectrum of hn contains discrete (bound) and continuum (scattering) states. 

First we obtain the bound states Cj, the solutions to eq. (0) with |a;| < m. The numerical 
method is sketched in Appendix A. We can use Levinson’s theorem to compute the number 
of bound states in each channel from the phase shifts. The phase shifts are computed from 
the S'-matrix, which in turn is extracted from solutions to second-order diherential equations 
obtained from the Dirac equation. Because we restrict our attention to backgrounds in the 
spherical ansatz, there are two conserved quantum numbers, grand spin and parity. The 
grand spin G is dehned as the vector sum of isospin and total angular momentum (orbital 
plus spin), and can be interpreted as a generalized angular momentum. The parity H is 
associated with space reflection in the usual way. 

We obtain second-order differential equations for the upper and lower components of 
the Dirac equation in the standard basis. After projecting onto a subspace with dehnite 
energy, grand spin and parity, we have two coupled second-order differential equations for 
two radial functions, gi{r) and ( 72 (^)- Together, the linearly independent solutions with 
incoming spherical waves in either of these channels dehne a two-channel scattering prob¬ 
lem. In the following we will suppress the labels cu, G, and H, which characterize this 
two-dimensional problem. The two linearly independent scattering boundary conditions are 
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labeled by {i,j} = 1,2 and are implemented as follows: At large r the solution has an 
outgoing wave ii i = j, and the radial wavefunction i ^ j vanishes. We summarize the two 
wavefunctions and two boundary conditions in matrix form, Qij{r) = g\^\r). We then write 
Q{r) as a multiplicative modification of the matrix solution to the free differential equations, 


Q{t) = 



= F{r)H{kr) 


where H is diagonal and can be expressed simply in terms of Hankel functions. 




k 


[X = 


G (^) 
0 


0 

h^G\x) 


and H_(x) = 


4+1 (a^) 
0 


hcliix 



(18) 


(19) 


for If = +(—1)*^ and If = —(—1)*^, respectively. For each value of grand spin, G, the parity 
quantum number dictates the values of i and f which enter H± = diag{hf\ In the 

channel with parity If = +(—)‘^ we have i = i' = G while for If = —(—1)*^ we have i = G + 1 
and i' = G — 1. Thus we have The elements of the 2x2 matrix F{r) satisfy second-order 
differential equations obtained from the Dirac equation. They are of the general form 


F” = — F' (1 + rL'(kr)) + -— 

r s Foj/m 


{F' + FL'ikr)) -VF + — [K, F] 


for upper and lower components respectively, where 


K = 


fiii + l) 0 
1 0 f (f -f 1) 


( 20 ) 


( 21 ) 


with £ and £' as above. V is the 2x2 matrix describing the coupling of the fermions to 
the Higgs background. The particular forms of V are listed in Appendix A. The matrix 
L = In H is the only remnant of the Hankel functions. 


L+{x) = 


In h 


( 1 ), 

G ' 


X) 


0 

In h^G\x) 


and L_(x) = 


In 


0 

In h^cl^ix) 


( 22 ) 


The elements of L'{x) = dL{x)/dx can be expressed as simple rational functions, which 
avoids any instability in the numerical treatment that would be caused by the oscillating 
Hankel functions. 

The 2x2 submatrix of the S-matrix can be constructed by superimposing solutions 
to eq. (^). First we normalize F by imposing the boundary conditions F{r ^ cx)) = 1 
and F'{r cx)) = 0. Given these boundary conditions, since the second-order differential 
equations for the gt are real, the scattering wavefunction can be written as 


Tsc = -F*{r)H*{kr) + F{r)H{kr)S (23) 

where S is the 2x2 submatrix of the F-matrix that we are seeking. Requiring that the 
scattering solution be regular at the origin yields 

S = \i-mH-^ikr)F-^ir)F*ir)H*ikr) . (24) 

r —>0 
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The quantity that enters the density of states is the total phase shift 


6{k) = — tr In S' = ^ lim tr In (^F ^(r)F*(r)^ 


( 25 ) 


from which H cancels because as r —> 0 the leading (singular) piece of H{kr) is real, i.e. 
lim^^o H*{kr)H~^(kr) = 1. The unitarity of S guarantees that equation (^) explicitly 
yields a real phase shift. 

Eq. (|25|) only gives the phase shift modulo vr. Of course, 5{k) should be a smooth function 
and vanish as fc —> oo. An efficient way to avoid spurious jumps by vr in the numerical 
calculation of 5{k) is to dehne 


(5(/c, r) = — tr In 

^ t 


■1-1, 


r)F*{r] 


(26) 


By construction 6{k) = lim^^o ?")• We then consider 


d6{k,r) 1 ^ 

dr 2i 


_ ^pp-^ 

dr dr 


(27) 


as an independent function to be included in the numerical routine that integrates the differ¬ 
ential equations for F, with the boundary condition lim^^oo ^(^, ^) = 0. Then hmr^o^(^)^) 
will then be a smooth function of k and go to zero as k ^ oo. 

We construct the Born series for F(r) in Appendix A. We introduce where 

ris and Up label the order in the expansion around s(r) = 1 and p{r) = 0 respectively. Then 
we hnd for the first two orders 


6^^\k) 

5^^\k) 


— lim tr 

2i r ->0 

— lim tr 

2i r ->0 


F<‘>"(r) - F6)(r) 


(28) 

(29) 


where 


^(1) = ^(1,0) ^ ^(0,1) pi2) ^ pi2fl) ^ ^(0,2) ^ ^(1,1) 

Subtracting these two terms from the full phase shift eliminates the quadratic divergence 
from the vacuum polarization energy. Eliminating the logarithmic divergence would be 
considerably more complicated because an expansion up to fourth order in n* -|- Up would be 
required.^ In Appendix B we introduce a simplified treatment for the logarithmic divergence. 
The resulting expression for the vacuum polarization energy, eq. (|^), then becomes 

®™c = -ii:(2G, + l)|£,|-|“^yPT^Ts(A.) + (31) 


^When restricting to field configurations with = v'^, two subtractions are sufficient 














m = i: (2G + 1) (iG,,,n(*:) - 4‘,l,n('=) “ Sal. 

G,cr,n 


n 


m 


k 


= 


- m 

^2 ^oo 


— arctan — + , „ „ , , 

m k k^ + / Jo 


J drr‘^ ^s(r)^ + p(r)^ — — 4(s(r) — 1)" 


(32) 


Tl^ Jo L J '■ 


Q j dx + x{l — x)q‘' 


In 


+ x{l — x)q^ 
m? — x{l — x)m 


H' 


m 

TT^ Jo 


dqq^p^{q)^m‘jj + 2711^ f da; 3 In f 1 — a;(l — x) 


m 


H 




— 2 In 1 + x(l 


--)4)]} 

J s J 


(33) 


( 4 ) _m.* (nf„ 


„ , „ + 6 / da; In 

OTT \ Jo 


1 — x{l — x) 


m 


H 




r/ \ 2 

xj drr"^ (^s(r)^ + p(r)^ — ij — 4(s(r) — 1)^ 


(34) 


where the last term in 5{k) implements the subtraction of the logarithmic divergence and 
is compensated by the terms in (see Appendix B). a denotes the sign of the energy 
eigenvalue, so that oj = a\/k^ + m?. We have introduced the Fourier transforms 


h{q) = 


' 2 sin gr /■“ , 2 

drr -('S(r) — 1) and p[q) = drr 

qr Jo 


cos qr sm qr 


qr 


(qr) 


p{r) . (35) 


Next we need to hnd the fermion number carried by the background held. We could 
consider an adiabatic transformation from the trivial background to the conhguration at 
hand, and count the levels that cross zero energy. In general, however, such a procedure 
is cumbersome. Instead, we apply the procedure developed in Ref. H. which is based 
on Levinson’s theorem. In each channel, we compare the number of positive energy bound 
states n^ch with the number of bound states that have left the positive continuum ^ ^G,+,n( 0 ). 
If one level originating in the positive energy continuum crosses zero, we will hnd that 
= “<^G,+,n( 0 ) — 1. In this case, the polarized vacuum carries charge {2G + 1) for each 
species. In general, the vacuum charge is given by 


G,n 


Qvac — “'^G,+,n(0) — rJ-Gu 

Ltt 


(+) 


--X!(2G' + 1) -,n(0) - 

G,n 


(36) 


The second equation rehects the equivalent counting procedure for negative energy states. 
We are interested in conhgurations with fermion number 1. If Qvac = 0 the fermion number is 
obtained by explicitly occupying a level. The lowest energy cost arises from hlling the level 
with the largest binding. If Qvac = 1 fhe polarized vacuum already provides the fermion 
number and none of the bound states needs to be explicitly occupied. 
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4. Numerical Analysis 

Our formalism is set up to allow consideration of an arbitrary background <f>(a;) of the 
form (^). However, as in all variational methods, we limit ourselves to variation of a few 
parameters in an ansatz motivated by physical considerations. We will scale energies and 
lengths in terms of the fermion mass m. We choose a four parameter ansatz 


s + ir ■ p = p(^) exp {ir ■ f 0 (^)) 

exp 


p (0 = 1 + ^1 

0(0 = -vr 


1 + 0 - 


w 


-1 


- 3 + 2e'=+/'^ 



(37) 


where ^ = mr, and the variational parameters are w, bi, 62 and c. Note that 0(0) = —vr, 
p( 0 ) = 1 + 61 , and both p — I and 0 go to zero exponentially as ^ > cx), since we expect a 

Yukawa tail. As long as p does not vanish, this background has winding number one. The 
width w is chosen such that 0(tc) = —vr/S. Furthermore, we have ensured that ^p(OI^=o “ 
0 , as the classical equations of motion require. 

In terms of p(0 and 6^(0, the classical energy eq. O is 


1 

m 


Eci(w,6i,62,c) 




(38) 


where pn = — • Then the total energy of the conhguration with fermion number 1 is 


—Etot 
m 


—£’ci(ta, bi, 62, c) + (1 - Qvac) ei{w, bi, 62, c) + Tvac(w, &i, &2, c) 


(39) 


where £’vac = E^a_c/rn, and ei = cui/m is the energy eigenvalue of the most strongly bound 
state. Note that for hxed pn, the coupling g appears only in the coefficient of the classical 
term. 

Conhgurations with 0(0) = —vr and 0(cxd) = 0 tend to strongly bind a state originating 
from the positive continuum in the = 0^ channel. For w large enough, this bound state 
will have crossed zero, causing the polarized vacuum charge to be Qvac = 1. In that case the 
level is not explicitly occupied and the corresponding term drops out of eq. (p9|). 


4.1 Sample Numerical Calculations 


For a given set of variational parameters, we hrst compute the phase shifts and perform the 
subtractions according to eq. which allows us to carry out the momentum integral in 


eq. (|3T[) . Using Levinson’s theorem, we then hnd the number of bound states in a given 
channel from 6G,a,u{k = 0 ), and use shooting to compute the bound state energies ej once 
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Figure 1: The total energy as a function of the Yukawa coupling constant g with mn = 0.35u. 


we know how many to look for. In terms of the scaled variables, the bound state and phase 
shift contributions in eq. (|3T|) depend on ansatz parameters, but not on model parameters. 
The dependence on model parameters is completely contained in Sd/g'^, /m and 

which are simple integrals involving the background helds. They can easily 
be obtained numerically for large regions in the model parameter space. Hence an efficient 
strategy is to choose a set of variational parameters and then consider the total energy as a 
function of the model parameters for that particular background conhguration. In Fig. |^we 
display a typical result for the total energy as a function of the Yukawa coupling g. 

The existence of conhgurations with total energy Etot/'m < 1 shows that there is a stable 
soliton whose energy is at most Etot- Apparently a sizable Yukawa coupling g is needed to 
obtain a stable soliton. However, as we will discuss later, our model is not reliable for such 
large Yukawa couplings because the Landau pole appears at an energy scale comparable to 
1/w. 

In Fig. I we display the total energy as a function of the depth parameters bi and 62 , 
for various values of the Yukawa coupling constant g and typical values of the remaining 
variational parameters w and c. We observe a shallow local minimum in the vicinity of 
bi = —0.8 for small and moderate values of g. However, at this minimum the total energy 
is larger than the mass m of the free fermion. For larger g, we obtain a total energy less 
than m for conhgurations with bi ^ 0. Conhgurations with 61 > 0 are more strongly bound, 
but the one fermion loop approximation fails for such conhgurations because the energy 
functional is not bounded from below (see Appendix C). We have therefore restricted the 
space of variational parameters to conhgurations for which the vacuum is stable to one loop. 

We observe a similar behavior for Etot as a function of 62 - There exists a local minimum for 
small and moderate g that does not yield a bound soliton. For large g, a bound soliton seems 
possible if 62 is big enough. In this case, the vacuum is stable at one loop for these values of 
the variational parameters. Note that when we hnd a marginally bound conhguration, the 
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Figure 2: Total energy as a function of the depth parameters bi and 62 with mjj = 0.35u 
and for various values of the Yukawa coupling constant g. 


Table 1: Comparison of the classical and renormalized vacuum polarization energy with the 
derivative expansion, cf. eq. (^Ol). 



w 

1.0 

2.0 

2.5 

3.0 

9 

= 5.0 

-0.127 

-0.024 

-0.008 

0.003 

9 

= 10.0 

-0.240 

-0.041 

-0.013 

0.005 

9 

= 15.0 

-0.287 

-0.048 

-0.015 

0.007 

9 

= 20.0 

-0.308 

-0.050 

-0.015 

0.006 


vacuum polarization contribution to the energy F^vac tends to almost exactly compensate for 
the gain from binding a single level. 


4.2 Comparison with the Derivative Expansion 


In order to check our computation of the vacuum polarization energy, in particular the 
simplihed treatment of the logarithmic divergence, we have compared our results with the 
derivative expansion. The relevant formulas for the derivative expansion are provided in Ap¬ 
pendix C. Denoting by Fgrad the energy computed to second order in the derivative expansion 
as obtained from eq. (|C. 6 |) , we display 


Ai = 


£c\ + Sv 


s, 


grad 


£cl + ^vac + 


grad 


(40) 


as a function of the width parameter w. The other variational parameters are kept constant 
at c = 1 . 0 , bi = —0.4 and 62 = 1-0. Also, we consider various values for the coupling constant 

9 - 

Where the derivative expansion is valid we expect Ai to go to zero. From Table ^ we 
conclude that our calculation agrees with the derivative expansion at large w, where the 
derivative expansion becomes exact. In this region, the derivative expansion is a good check 
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Table 2: Comparison of the total energy with the gradient expansion, cf. eq. (^T]) . 


w 

1.0 

2.0 

2.5 

3.0 

9= 5.0 

0.168 

0.034 

0.018 

0.012 

g = 10.0 

0.257 

0.056 

0.029 

0.019 

g = 15.0 

0.285 

0.063 

0.033 

0.022 

g = 20.0 

0.297 

0.066 

0.034 

0.023 


on our computation of the vacuum polarization energy. In particular, it checks our handling 
of renormalization, since renormalization effects are included in the derivative expansion 
in the standard way. On the other hand, it is clear from Table |I] that the second-order 
derivative expansion cannot be trusted for w ~ 1.0, i.e. for background configurations 
whose extension is close to the Compton wavelength of the fermion. One might imagine 
improving the derivative expansion by including the effect of an explicitly occupied level. 


^2 




tot 


s, 


grad 


£tot + 


grad 


(41) 


because the derivative expansion to Qvac, the topological charge, suggests that the back¬ 
ground configuration carries fermion number regardless of the value for w. However, this 
modification does not give any better agreement at small w, as can be observed from Ta¬ 
ble In general, the inclusion of the explicitly occupied level tends to change the sign of 
the relative deviation. 


4.3 The Landau Pole 

From these studies one might conclude that a soliton takes over the role of the lightest 
fermion once the Yukawa coupling constant g becomes large enough. At large g, the positive 
contribution to the total energy from Td in eq. (p8D, which disfavors the soliton, decreases 
quickly for large g. However, for large couplings the model itself becomes ill-defined. Since 
the model is not asymptotically free, it has a Landau singularity in the ultraviolet, reflecting 
new dynamics at some cutoff scale. Thus the Landau pole sets a minimum distance scale 
below which the model is not consistent. Solitons that are large compared to this scale are 
relatively insensitive to the unknown dynamics at the cutoff scale, but solitons whose size is 
comparable to this scale cannot be trusted. In this section we will discuss the emergence of 
the Landau pole and estimate its effect on the vacuum polarization energy by comparing the 
present results with a calculation that removes this pole. Although this removal is somewhat 
ad hoc, it nevertheless provides some insight into the reliability of the computations in case 
of large g. 

Denoting the Fourier transforms of s(r) — 1 and p(r) by h{q) and p(^, respectively, the 
contribution of the two-point function to the total energy can be written as 

^2 = y y ■ p{-^] (42) 
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where 


G-\q^) 


V 

y 


(g2 - ml) [ 1 + 


47r^ 


47J-2 


6 dx 


rn? — a;(l — x)q^ 


In 


IT 9^1 + 


9 


+ 


9 


47r^ / 47r^ 


Q dx 
Jo 


m? — a;(l — x)q^ 
m? — x(l — x)ml 

— x{l — x)q^ 


(43) 


In 


rn? — x(l — x)q^ 
T 

H 


im? — x(l — x)m 


+2m / dx 
Jo 


3 In 1 — x(l — x) 


m 


H 




— 2 In 1 — x{l — x)- 




(44) 


which includes classical, loop, and counterterm contributions. Gp{ql has a pole (the Landau 
ghost pole) at space-like momentum —ml with residue Zq. The pole location is easily 
obtained numerically from the condition G~^{q‘^ = —ml) = 0. In the vicinity of q^ ~ ~^l 
we have the expansion 


with 


Gp^{q^) = ^ (f + ml) +0(f + m?G)‘ 


1 


'^Gp^i^ql 


V 

y 


1 - 


yl 


O 9 

q^=—mQ 

^ N, -I-a;(l — a;)m^ 

6 / dxx{l — x) In ^ 

Jo 


—4 / dx 


m^ — a;(l 
x{l — x) 


x]m 


H 


0 -|- x(l — x)ml 


(45) 


(46) 


The existence of this pole yields an unphysical negative contribution to the total energy at 
large space-like momenta, or equivalently for narrow background held conhgurations. Based 
on the Kallen-Lehmann representation for the two-point function, the authors of Ref. [ 1^ 
suggested a procedure to eliminate the Landau pole while maintaining chiral symmetry. 
They replace eq. (^2]) with 


Eo = 




(47) 


where Ap(g^) = Gp{ql — Zq! (g^ + ml) removes the Landau pole. We can easily adopt this 
procedure since we have already extracted the loop and counterterm contributions from the 
two-point function in eq. ([531). That is, we replace + E^l + E^^l by E 2 + T'dwith 


-E. 

m 


( 3 , 4 ) 


7^ 

2g2 



2 (s - 1) (s - 1)^ +p‘^ 


4(.-1)2 



(48) 


In Fig. ^ we show the effect of this replacement as a function of the width parameter w for 
g = 12, which is in the region where a bound soliton can occur. For small w, we observe that 
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Figure 3: The Landau pole subtraction, for run = 0.35u and g = 12. 


the original computation, eq. (|^, gives a negative contribution. However, for small w there 
are only weakly bound states and the vacuum is almost undistorted. The classical energy is 
also small since g is large. Hence the total energy is dominated by the renormalized Feynman 
diagram contribution which can be negative due to the Landau pole. Thus for 

small w and large g, the Landau pole dominates the binding of the soliton, and, even worse, 
the total energy could be negative, reflecting an unphysical vacuum instability |^. Using 
the above prescription to eliminate this pole, the total energy turns out to be positive for 
all values of w, so the instability is removed. As can be seen from Fig. ^ for sensible w 
this prescription increases the total energy by about 0.25m for the parameters chosen, which 
unbinds the soliton. We conclude that the solitons found at large g are principally bound 
by unphysical effects associated with the Landau singularity, and not by reliable dynamical 
properties of the model. 


4.4 Scalar Backgrounds 


We conclude the numerical analysis with a calculation of the total energy when only a scalar 
background held s(r) is present. Our goal is a brief comparison with the results of Ref. p!7[] , 
rather than a complete study. For p = 0 the Dirac Hamiltonian is charge conjugation 
invariant. Hence the charge carried by the background held is zero and we must explicitly 
occupy the most strongly bound state. In Fig. § we show typical results of the numerical 
calculation for the total energy as a function of the coupling constant. The hgure shows that 
a slightly bound soliton emerges even for modest values of the Yukawa coupling. Its energy 
is up to 5% less than that of a fermion in the background of the translationally invariant 
scalar held. In this case, the Landau ghost singularity in Gp (c/. eq. (1^)) is irrelevant and 
we may trust the calculation even for small w. Furthermore, the second-order derivative 
expansion deviates from the full calculation by only a fraction of a percent even at moderate 
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Figure 4: The total energy with only the scalar background as a function of the Yukawa 
coupling constant g with rriH = 0.35u. Note that the variational parameter c is irrelevant 
without pseudoscaiar heids. 


values of the width parameter, w ~ 1.5. Plotting the results from the derivative expansion 
would yield indistinguishable curves in Fig. Our approach thus conhrms the hndings of 
Ref. |]T3, which used the derivative expansion to hud a slightly bound soliton. 


5. Summary and Outlook 

We have performed a quantitative search for Higgs solitons in a theory with chiral 
fermions. In the analogous model in one spatial dimension, such solitons exist and are 
strongly bound for a wide range of coupling constants. In three spatial dimensions, however, 
we did not hnd any region where the soliton binding was strong enough to be convincing. 
Twisted Higgs background heids do strongly bind a fermion level, but it is necessary to add 
the renormalized energy due to the distortion of the fermion vacuum. We have developed 
new methods to regularize, renormalize, and compute the corresponding vacuum polarization 
energy. In this model, we hnd that the vacuum polarization tends to cancel the contribu¬ 
tion of the strongly bound fermion level. The total one-loop energy overcomes the classical 
energy only for large Yukawa couplings, where the theory has an unphysical Landau pole. 
Hence energetically favored Higgs solitons observed for large Yukawa couplings cannot be 
interpreted as reliable predictions of the model. 

Nevertheless, these hndings do not rule out the existence of Higgs solitons within the 
Standard Model. In particular, the inclusion of the gauge heids could be critical to soliton 
formation. This point of view is motivated by a number of considerations: 

• Expanding the variational ansatz to include gauge heids can only decrease the energy 
of a conhguration. In this respect, the 3-1-1 dimensional model is diherent from the 
1-1-1 dimensional model where we did hnd a soliton, because in one dimension gauge 
heids do not add any new interactions. 
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• Gauge fields are essential to the anomaly arguments underlying the decoupling results. 

• In the gauge theory, there is a sphaleron barrier with height 2Tiv/gw, where v is the 
Higgs vacuum expectation value and gw is the gauge coupling. If the fermion has 
a Yukawa coupling gp such that its perturbative mass gpv is much larger than the 
sphaleron height, it has an unsuppressed decay mode over the sphaleron. Thus the 
ordinary fermion states are unstable and are nowhere to be found in the spectrum of 
the fermion Hamiltonian. The creation of a soliton state with mass below the sphaleron 
energy would allow the theory to remain gauge invariant after decoupling as required 
byi. 

Work is now underway to generalize this calculation to include the gauge helds. 
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Appendix A: Dirac equations 


In this Appendix we present the Erst- and second-order Dirac equations used in Section 3 
for fermions in a static background held in the spherical ansatz, eq. ( 0 )- In order to solve the 
Dirac equation, eqs. (BO). we begin with spinors that are eigenstates of parity and total 
grand spin [^, where grand spin G is the sum of total angular momentum j = I + and 
isospin ir. For a given grand spin quantum number G with ; 2 -component M, we will End 
the bound state and scattering wavefunctions in terms of the spherical harmonic functions 
yj/{x) with j = G ± I and £ = j ± These are two-component spinors in both spin and 
isospin space. In the following we will suppress the label M for the z-component. While 
grand spin is conserved, the background field in eq. (|^ mixes states with difierent total 
angular momentum j. For the channels with parity H = -l-(—the spinor that diagonalizes 
eq. (pRl) reads 


_ ( ^9i{r)yG+id^) \.( ig2{r)yG-i^,Gi£) \ 

^ [h{r)yG^,G^i{£)) ^ [-f2{r)yG-lG-i{£) J 

The spinor with opposite parity, —is parameterized as 




(f) 


( m{r)yG+^,G+iix)\ ( 

V -fi{r)yG+i,Gi^) J [ f2(r)yG-lG(£) J 


(A.l) 


(A.2) 
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Note that for convenience we have again suppressed both the grand spin and parity labels 
for the radial functions Qi and /j. 

The matrix elements of the operator r • i. between various yji can be found in the 
literature Then the Dirac equation can be expressed as coupled hrst-order differential 
equations for the radial functions Qi and ft. In the 11 = +(—)*^ parity channel, we have 


0 = — —gi + (ms + u!)fi — 2Q^ I (^1 ~ 2 y ^ G{G + l)g2^ 

0 = <72 “I- - —fi's ~ + <^)/2 + 2G + 1 

0 = /i H- -—fi + {ms — u!)gi + 2Q^ i 2^ G{G + l)/2^ 

0 = /2- - —/2 — {ms — u!)g2 — 2G + I ~ 2^ G{G + l)/i^ (A.3) 

where a prime indicates a derivative with respect to r. In the 11 = —(—)*^ parity channel, 
we have 

0 = S'! H- -—gi — {ms + uj)fi — 2 q^i ~ ‘^\fG{G~^^g‘^ 

0 = ^/2- - — g2 + {ms + oj)f2 + 2 G ^1 

0 = /i — —fi — {ms — uj)gi + 2^ ^ ^ G{G + l)/2^ 

0 =/2 H- - —/2 + {ms — uj)g2 ~ 2 G +I ~ ^ l)/i^ • (^- 4 ) 


We can use these differential equations to obtain the bound state solutions with |a;| < m 
using ordinary shooting methods. For a number of cases, we have verihed the solutions 
numerically by diagonalizing eq. (pTl) in a spherical cavity As discussed in Section 3, we 


also require the second-order equations for either the upper {gi) or lower {fi) components of 
the Dirac spinor. These components will be collected into the matrix F{r)H{kr) dehned in 
eq. ([I8|). In Section 3 we chose to work with the upper components only. Here we will discuss 
both upper and lower components, Fu{k) and Ffk). Substituting this parameterization into 
the second-order equations for the upper or lower components of the Dirac equation and 
multiplying by the inverse of H{kr) from the right, we obtain for the upper components in 
the channel with H = -f-(—)‘^ 


0 = F: + -Fl\l + rL'4kr) 


s + ^ L 




(A.5) 


with F = Fu entering eq. (^0|). The matrix L^{kr) is dehned in eq. (|^) . For later conve¬ 
nience we have added the fermion mass m as an argument of the potential matrix. 


Vn{r, m) = —rnf s^ + — 1 


+ 


G 


r s + - 

m 


mp' 
2G + 1 


2{G + l)mp 
r 2G + 1 ^ 


mp s' 

2GTTsT^ 

m 
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V 22 {r, m) = —m^ — 1 


G + 1 s' 


+ 


mp' 2 Gmp 


mp 


r s + ^ 2G + 1 r2G + l 2G + 1 s + ^ 


Vi2{r,m) = V2i{r,m) = 


2^G{G + l) 
2G + 1 


mp + mp - 


r s + 


(A.6) 


Similarly, we find the second-order eqnation for the lower components in the channel with 
n = +(-)G is 


0 = F!' + -F! 1 -I- rL'_(kr) 

r L 


where 


hhii (r, m) = —m^ — 1 


W 22 {r,m) =—m s +p — 1 


s- ^ L 


G + 2 s' 


F; + FiL'_{kr) + W{r, m)Fi - [K, Fi] (A.7) 


+ 


mp' 2{G + l)mp 


mp 


r s-^ 2G + 1 r 2G + 1 2G' + ls-^ 


+ 


G-I s' 


mp' 2 Gmp ^ mp s' 


r s-^ 2G + 1 r2G + l 2G + 1 s - ^ 

m m 


Wi2{r,m) = W2i{r,m) = 


2^G{G + 1) 

2GTi 


/l s' 
mp' — mp —h 


r s + 


(A.8) 


In this case we have two different orbital angnlar momenta, leading to the commutator term 
K = (l/r^) diag {(G + 1)(G + 2), G{G — 1)}. Next we write down the differential equations 
in the channel with 11 = ■ Using dehnitions analogous to eqs. (|l^) we hnd for the 

upper components 


0 = F" + -F,^fl+rL'_(A;r) 




F' + FuL'_{kr)] + W{r, -m)F^ - [F, F„] (A.9) 


while the lower components obey 


0 = F” + -F/ [l + rL'^{kr) 


s -L 

m 


F/ + FiL'^{kr) \ + U(r, -m)Fi . 


(A.IO) 


We must study the Born series in order to obtain the subtraction terms in eq. (^). We 
expand around the free solution, 

Fu(r) = 1 + + 9,F<“’'>(r) + 9!q"'“V) + + 3-3pFF\r) ... (A.ll) 

and similarly for F;. The expansion parameters Pp and Qs are dehned by 

<I)(x) = n [1 -I- gsis{r) — 1) -I- ippT ■ xp{r)] (A.12) 


where we have expanded the equations of motion and the potential matrices V and W in 
terms of the artihcial coupling constants Qs and Pp. Having obtained the second-order Dirac 
equations at the desired order in the couplings, we also expand the dehning equation for the 
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phase shifts, eq. (^), in these constants. Then, for example, and will contribute 

to (^G.p.TT- However, we observe the relations 

vn, 2 i+i^^^ m) <—^ -m) and m) <—^ -m) (A. 13) 


when u <—»• —oj. As a result, once we sum over parity channels, the positive and negative 
energy pieces in the vacuum polarization energy calculation cancel for all odd powers of the 
pseudoscalar field. Of course, this cancellation just reflects parity invariance. 

For scattering solutions, we have |a;| > m and hence the second-order equations are 
nonsingular as long as |s| < 1. As explained in the main part of the paper we integrate these 
second-order equations from r—>'Ootor = 0. Atr—>cx)we have s = 1 and commonly s 
changes sign at some intermediate point, say tq. For cn > 0 it is hence appropriate to use 
the second-order differential equations for the upper components, which will be the larger 
ones. Eventually, however, s may become less than minus one. In that case we will have 
singularities when using the upper components. Then it would be more appropriate to 
employ the lower components as these will be the larger ones for r < tq. For a; < 0 the 
situation is reversed. We therefore switch between these two components at r = tq using 
the hrst-order Dirac equations. In the parity (—channel the switch from upper to lower 
components is given by 

Fi = (f: + F„L'+ + MFu) H+HZ^ (A. 14) 

uj -I- ms ^ ^ 

F; = -FiL'_ - MFi + (w - ms)a:iF^H+HZ^ • (A. 15) 


In this case the upper components are obtained by integrating from inhnity to tq and Fi in 
eq. (A-15|) is given by eq. (A-14|) . The matrices M and M are defined by 


mp 


I GlG+l) 


M = 


r _^+1 

V/G(G+1) 

^ 2G+1 ''''F 


2G-I-1 


mp 


G+2 


and M = 


+ 


mp 

2G+1 


G+1 I 
r 2G+1 


mp 


, \/G(G+1) 
' 2G-I-1 


mp 


I ^/G{G+T) 

‘ 2G+1 

G—1 _ mp 


mp 


2G+1 


(A.16) 

where the Hankel function matrices H± are dehned in eq. (R). The switch from lower to 
upper components is given by 

as [FI + FuLZ + 


F, = 


UJ — ms 

K = -FuL'^ 


(f/ + FuL'_ + MFi) H.H-^ 


In the parity 


MFi — {uj + ms)asFuH_H_^_^ . 

-)‘^ channels the switch from upper to lower components is given by 


(A.17) 

(A.18) 


1 


-as (f: + F^L'_ + NF^) 


UJ + ms 

F/ = -FiL'^ - NFi -{uj- ms)asFuH_H-^ 


(A.19) 

(A.20) 


with 


N = 


G+2 

r 


mp 

2GTT 


I VG(G-I-I) 
‘ 2G-I-1 


mp 


I a/G(G- l-l) 
‘ 2G-I-1 

-^ + 


mp 


and N = 


G I mp 
r ' 2G+1 


I a/G(G+1) 


mp 

2G+T 


I \/ G(G+1) 
‘ 2G-I-1 


mp 


- 2G+1 -^P 

G-l-1 _ mp 
r 2G-I-1 


(A.21) 
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and the transformation from lower to npper components is 


Fu = - (f; + F^L'^ + NFi) H+HZ^ (A.22) 

oj — ms ^ ^ 

Fi = -F^L'_ - NFi + (a; + ms)a^F^H+HZ^ . (A.23) 


Finally, we note that the reqnired ratios of Hankel fnnctions can also be expressed as rational 
fnnctions 


H+HZ^ 

H_H-^ 


a^L'_ H— 
r 

—a^L'_ + 


G + 2 0 \ 

0 G-l) 

G 0 Y 
0 G + l)_ 


(A.24) 


Appendix B: Fake Boson Field 

In this Appendix, we describe the special treatment of the logarithmic divergence of 
Feynman diagrams. In particnlar we will explain and provide nnmerical evidence for the nse 
of the simplihed form in eqs. (|3^) and (0). 


1. Discussion 

In a fermion loop calcnlation, diagrams with one or two external lines are qnadratically 
divergent and those with three or fonr are logarithmically divergent. Hence we will have to 
take A^max = 4 in eq. Althongh we can straightforwardly compnte the corresponding 
Born terms as in eqs. (^81) and (^), the eqnivalent Feynman diagrams are difficnlt to com¬ 
pnte nnmericallyQ when the external helds are coordinate-dependent. All we really need to 
do, however, is to regnlate the momentnm integral in eq. by snbtracting an appropriate 
expression from the integrand and ensnre that we add back in exactly the same qnantity. 
The latter qnantity shonld have a divergent piece that can easily be canceled by the conn- 
terterms. Boson loops have a mnch simpler divergence strnctnre: the logarithmic divergence 
corresponds to a Feynman diagram that is only second order in the external lines. Higher- 
order boson loop diagrams are hnite. We can therefore simplify the regnlarization of the 
logarithmic divergence of the fermion loop signihcantly by snbtracting and adding back the 
contribntions of an eqnivalent boson. This boson is completely artihcial to the model so 
we will call it the “fake boson held.” We impose the condition that at one loop the fake 
boson model generates the same regnlated logarithmic divergence as the original fermion 
loop. We then snbtract the associated second-order Born phase shift from the fermion phase 
shifts and add it back in as a regnlated Feynman diagram. By constrnction, its divergent 
piece is canceled by the connterterm contribntions in eq. This fake boson method does 

^The fourth-order diagram requires a nine-dimension integral, not including Fourier transforming the 
background fields. 
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not give the full Feynman diagrams of the fermion loop, so the approach is not suitable to 
determine counterterm coefficients in a specihed renormalization prescription. It can only 
be used once these coefficients are known. In our model we have uniquely determined the 
counterterms from the hrst- and second-order terms in the expansion of the fermion loop, 
which we computed exactly. We can then apply the fake boson method to the third- and 
fourth-order terms, which would otherwise be difficult to evaluate. We can extract the local 
piece of a Feynman diagram by setting the external momenta to zero. An expansion in the 
external momenta then shows that for a second-order boson diagram only this local piece 
diverges. We will identify the local piece in the second-order boson diagram as a “limiting 
function” to the the second-order Born approximation to the phase shift. This procedure 
provides the simplihed expression used in eqs. ( 0 )-®. 


2. Equivalent Boson 

We begin by considering the second-order Born approximation to a boson loop. We 
consider the spherically symmetric problem 


dr'^ 


ue{r) + 


'+ 1 ) 


gV{r) 


ue{r) = k‘^ue{r) 


(B.l) 


discussed in Ref. [^. The coupling constant is a bookkeeping device which at the end 
we will take to be unity. Solving for the complex function j3^{k,r) in the ansatz ug{r) = 
exp{2i(3£{k,r))rhf\kr) with the boundary conditions (3£{k,oo) = (3£{k,oo) = 0 yields the 
phase shifts [|T^ 

6e{k) =-2Re(3e{k,0) . (B.2) 

The differential equation for (3£{k,r) is non-linear and the expansion (3£{k,r) = g(3f\k,r) + 
g‘^(5f\k,r) + ... yields the various orders of the Born series to the phase shifts = 
—2Re/3£"'^(fc, 0) by iteratively solving the differential equations for (3^"'\k,r). V{r) provides 
the source term for f3^^\k,r). Formally, the quadratically and logarithmically divergent 
contributions to the vacuum polarization energy are contained in 


i:(2« + 

Jo Ti dk ^ 

JQ 7T dn^ £ 

Now let us consider two potentials, Vi(r) and V 2 {r), which are related by 

POO POO 

/ drr‘^V^{r) = / drr'^V^ir) 

Jo Jo 


(B.3) 


(B.4) 


where we also allow for different masses mi and m 2 of the boson held. The dispersion relation 
uj = \Jk‘^ + mf is the only place where a dependence on the mass appears since eq. ( p.l|) 
does not contain the mass parameter explicitly. 
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Table B.l: Comparison of second-order contributions to the vacuum polarization energy. 


d 

w 

m2/mi 

Ay 

Ap 

-2.0 

2.0 

1.0 

0.088 

0.087 

-2.0 

2.0 

0.5 

6.915 

6.916 

-2.0 

2.0 

2.0 

-7.309 

-7.312 

2.0 

0.8 

1.0 

0.010 

0.010 

2.0 

1.5 

1.0 

0.023 

0.023 

3.0 

2.0 

1.0 

0.011 

0.011 


Since the logarithmically divergent counterterm only depends on the potential through 
the quantity / d^xV{x)‘^, it will be identical for both potentials. Therefore the difference of 
the second-order Feynman diagrams 


Ap — 


dq 


(dvr) 2 Jo (27r)2 

1 /■°° dqq^ 

(47r)2 Jo (27r)^ 


Vi{q)V,{-q) 


-2 + 


\ 




q- 


1 + ^-1 
9 


V2{q)V2{-q) 


-2 + 


\ 


1 + 




In 


l + ipi + l 


^ 4m2 _ ^ 


(B.6) 


will be finite. Here Vi{q) denotes the Fourier transform of Vi{r). Since we can identify orders 
in the Feynman diagrams with those in the Born series for the vacuum polarization energy, 
Ap should be identical to 


Ay = E, 


( 2 ) 

caSjl 


E. 


( 2 ) 

cas,2 


(B.6) 


where the difference is to be taken under the integral in the second equation of ( |B.3|) . The 
second subscript, i = 1,2, refers to the potential Vi. In Table |B.1| we compare Ap and Ay 
for two Gaussian-type choices 


V^{r) = + 2de-’''/2-^ and V2{r) = 


(B.7) 


Here d and w are variational parameters and the coefficient C is hxed by eq. ( [B.d|) . We 
observe that the differences Ay and Ap agree within the numerical accuracy even though 
either of them may be quite large in magnitude, especially when the two masses are taken 
to be different. We conclude that we can subtract the second Born approximation and add 
back in the second-order Feynman diagram associated with bosonic fluctuations about Vi to 
regulate the logarithmic divergences encountered in the study of any other problem with the 
same divergences. The method we employ in the Yukawa model is the generalization of this 
procedure to the case of fermions. Note that by considering the second-order Dirac equation 
(c/. Appendix A) the phase shift calculation in the Yukawa model is essentially that of a 
boson held with derivative interactions. 
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3. Limiting Function 


Next we extract a local contribution containing the logarithmic divergence. We will 
manipulate this expression so that it can be substituted into the phase shift formula for the 
vacuum polarization energy. This procedure leads to further simplihcations for evaluating 
the fermion vacuum polarization energy. Since these manipulations involve divergent objects 
they are not rigorous results. However, we will be able to verify their validity numerically 
for a specihc background potentials. This check is sufficient to justify the use of these 
simplifications in the Yukawa model because we can always revert to that specific potential 
using the arguments of the previous subsection. 

We formally identify the local contribution by setting the external momenta to zero in 
the second-order Feynman diagram 




1 m 

— arctan — -|- 
m k 


W| 

^2 _|_ ^2 J 


(B.8) 


where we have carried out finite and angular integrals. Nevertheless, these manipulations 
are formal since they involve the logarithmically divergent fc-integral. However, so far we 
have only manipulated the integrand. It is worthwhile to note that the k dependent function 
in the last integral equals that of the last term in eq. (^). We subtract the local contribution 
from the full second-order Feynman diagram, giving 


= 


1 r°° dqq^ 

[AttY Jo (27r)2 


V{q)V{-q) 




1 + ^+1 


1 + ^ 


(B.9) 


which, of course, is finite for the potentials of the type listed in eq. (|B.7|) . Similarly we can 
define a finite second-order energy by subtracting the formal expression from the second 
Born approximant to the vacuum polarization energy 


= 


dk 




d 




/o vr dk 

where the limiting function is given by 




(B.IO) 


SiS,. tct.(k) = -o / drr^V^(r) i —arctan^ + 
8 JO \m k 


k 


k2 _|_ ^2 


(B.ll) 


In Table [B.2| we compare AEp'^ and AFy^ for the background potential Vi{r) given in 
eq. O- The particle in the loop of the local contribution does not need to have the same 
mass as that in the full Feynman diagram. Within the numerical accuracy we do not find 
any differences. 
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Table B.2: Comparison of the second-order Feynman diagram including the local subtraction 
with the corresponding expression for the vacuum polarization energy. 


d 

w 

m 2 /mi 



2.0 

0.8 

1.0 

0.001 

0.001 

2.0 

1.5 

1.0 

0.002 

0.002 

3.0 

2.0 

2.0 

0.113 

0.112 

1.0 

1.5 

0.8 

0.069 

0.069 

1.0 

1.5 

1.2 

-0.036 

-0.036 


Appendix C: Derivative expansion 


Using the techniques developed in Ref. we compute the two leading orders of the 
derivative expansion for the fermion determinant 


r(<F) = —i Tr In {i^ — . 


(C.l) 


First we compute the effective potential. In dimensional regularization, 


r('i>) 


2 

(47r)2 



47rp,2 




(C.2) 


where /i is the scale introduced to render g dimensionless in d = 4 — 2e dimensions. In order 
to extract the contribution with two derivatives acting on <F we parameterize $ = <Fo -|- 
with <Fo constant. Since 

F(<ho + (5$) = -U(<ho + d<h) + ^Fg](<Fo)a^d<Fa^d<F + 0{d5^f (C.3) 

it is sufficient to expand F(<Fo -|- 5$) to quadratic order in both <5$ and the momentum of 
its Fourier transformation. Returning to conhguration space this yields 


i 

6 




(g2 - 




(C.4) 


Now we replace <Fo by <F and 9^5$ by and treat the loop integral in dimensional regu¬ 
larization, 



(7 - 7 - ^ - In tr^ 

+^tr^ j 


(C.5) 
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Combining the expressions eq. ( |C.2|) and (|C.5|) with the counterterms computed in Section 
3 yields the hnal result for the derivative expansion up to second order, 


£(2)^ 


(dvr) 


C 2 - In 


l<hl 


triT 


+^trF j I 




IbTT^ 


i$r-u^-2i<hrin 


l<hl 


where 


mir 3 /•! , , 

"» = 4^?+2/„ 

C 2 = 6 / dxx{l — x) In 
Jo 


1 — x{l — x) 


m 


H 




1 — x{l — x) 


m 


H 




(C.6) 


We observe that for conhgurations with the sum of the classical energy and 

the contribution computed from eq. (|C.6|) can become negative. Such conhgurations will 
destabilize the vacuum. Numerically, we have verihed that this behavior also emerges in 
the full calculation for the vacuum polarization energy, which goes beyond the derivative 
expansion. Certainly this is an artifact of the one-loop approximation, implying that within 
this approximation we may not consider such conhgurations. 
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